Effects of dynamical phases in Shor's factoring algorithm with operational delays 
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Ideal quantum algorithms usually assume that quantum computing is performed continuously by a sequence 
of unitary transformations. However, there always exist idle finite time intervals between consecutive operations 
in a realistic quantum computing process. During these delays, coherent "errors" will accumulate from the 
dynamical phases of the superposed wave functions. Here we explore the sensitivity of Shor's quantum factoring 
algorithm to such errors. Our results clearly show a severe sensitivity of Shor's factorization algorithm to 
the presence of delay times between successive unitary transformations. Specifically, in the presence of these 
coherent "errors", the probability of obtaining the correct answer decreases exponentially with the number of 
qubits of the work register. A particularly simple phase-matching approach is proposed in this paper to avoid or 
suppress these coherent errors when using Shor's algorithm to factorize integers. The robustness of this phase- 
matching condition is evaluated analytically or numerically for the factorization of several integers: 4, 15, 21, 
and 33. 



PACS numbers: 03.67 .Lx 

I. INTRODUCTION 

Building a practical quantum information processor has at- 
tracted considerable interest during the past decade 1 1]. With 
the resources provided by quantum mechanics, such as su- 
perposition and entanglement, a quantum computer could 
achieve a significant speedup for certain computational tasks. 
The most prominent example is Shor's factoring algorithm 
Ql^l, which allows an exponential speedup over the known 
classical algorithms. The proposed quantum algorithms are 
constructed assuming that all quantum operations can be per- 
formed precisely. In reality, any physical realization of such a 
computing process must treat various errors arising from var- 
ious noise and imperfections (see, e.g., 01). Physically, these 
errors can be distinguished into two different kinds: incoher- 
ent and coherent errors. The incoherent errors originate from 
the coupling of the quantum information processor to an un- 
controllable external environment, which is stochastic, and re- 
sults in decoherence. Coherent errors usually arise from non- 
ideal quantum gates which lead to unitary but non-ideal tem- 
poral evolutions of a quantum computer. So far, most previ- 
ous works (see, e.g., 0111 Gal) have been concerned with 
quantum errors arising from the decoherence due to interac- 
tions with the external environment and external operational 
imperfections. Here, we focus instead on internal ones. The 
coherent errors we consider here are related to the intrinsic 
dynamical evolution of the qubits between operations. 

A quantum computing process generally consists of a se- 
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quence of quantum unitary operations. These transformations 
are usually applied to the superposition states so that the quan- 
tum computer evolves from an input initial state to the de- 
sired final state. If the two qubit levels have different ener- 
gies, as it is usually the case, the superposition wave function 
of the quantum register undergoes fast coherent oscillations 
during the finite time delay between two consecutive opera- 
tions. These oscillation, if not controlled, can spoil the correct 
computational results expected from the ideal quantum algo- 
rithms, where operational delays are neglected. 

In principle, these coherent errors can be either (1) avoided 
by tuning the relevant energy splittings of the qubits to zero 
II (I Hill : or (2) eliminated by introducing a "natural" phase 
induced by using a stable continuous reference oscillation for 
each quantum transition in the computing process ifll . Ex- 
perimentally, for example in NMR systems (see, e.g., O lOO . 
these errors were usually corrected by introducing two addi- 
tional operations before and after the delay to reverse each 
undesired free evolution. 

In this paper we perform a quantitative assessment of the 
effects of the dynamical phases in Shor's algorithm by real- 
istically assuming that operational delays, between succes- 
sive unitary transformations, exist throughout the computa- 
tion. We explore a phase-matching approach to deal with 
the dynamical phase problem. We show that coherent "er- 
rors" due to these phases, acquired by the dynamical evolu- 
tion of the superposed wavefunction during the operational 
delays, may be avoided by properly setting the total delay. 
We then carefully evaluate the robustness of such a phase- 
matching condition, focusing on its dependence on the num- 
ber of qubits, the length of the delay, and the fluctuations in the 
qubit energy splitting. Our discussions are in the context of 



Shor's algorithm, but can be extended to other quantum algo- 
rithms, such as the phase estimation and other algorithms 11411 . 
For simplicity and clarity, here we assume that the influence 
of the environmental decoherence and the gate imperfections 
on the computing process are negligible. 

The paper is organized as follows. In Section II, we present 
a decomposition of Shor's algorithm and explain how we in- 
corporate the dynamical phases into the realization of this al- 
gorithm. The usual decompositions of quantum algorithms 
into consecutive elementary gates are strictly limited by the 
short decoherence time. Here, we reconstruct the standard 
Shor's algorithm out of four functional unitary transforma- 
tions, and only consider the operational delays between these 
larger building blocks. We assume that each block can be ex- 
actly performed by only one-time evolution as a multi-qubit 
gate (see, e.g., 1151 lion , avoiding the existing idle time in- 
side it. It is shown that the effects of dynamical phases are not 
negligible, even in this primary or "coarse-grained" decompo- 
sition. In Section III, we numerically evaluate several exam- 
ples to illustrate the phase-matching condition, and establish a 
clear relationship between this condition and the equivalence 
of the Schrodinger picture and the interaction picture descrip- 
tion of a physical system. We also demonstrate the robust- 
ness of the phase-matching condition by varying the number 
of qubits involved, the delay duration, and distribution of qubit 
energy splitting. Finally, in section IV we present some con- 
clusions and discussions from our numerical studies. 



II. FOUR-BLOCK DECOMPOSITION OF SHOR'S 
ALGORITHM WITH OPERATIONAL DELAYS 
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FIG. 1: Quantum circuit for implementing Shor's algorithm with 
time delays t+ ( j = 1,2,3) between the successive operations. 
Here H refers to a Hadamard gate, while F refers to a quantum 
Fourier transformation. Each block operation is assumed to be ex- 
actly performed in a very short time interval e (so that phases accu- 
mulated during the operations are either accounted for by the opera- 
tions themselves or simply neglected). 



computational initial state of the system becomes: 

9-1 



Here, the subindex W stands for the work register state, and 
the subindex A for the auxiliary register. After a finite time 
delay n, and right before the second unitary transformation is 
applied, the initial state |W(0)) of the whole system evolves to 



We study the dynamical phase problem in the context of 
Shor's factoring algorithm. In Shor's algorithm |2], the fac- 
torization of a given number N is based on calculating the 
period of the function f(x) = a x modiV quantum mechani- 
cally for a randomly selected number a (1 < a < N) coprime 
with N. Here y mod N is the remainder when y is divided 
by N, The order r of a mod N is the smallest integer r such 
that a r modiV = 1. Once r is known, factors of N are ob- 
tained by calculating the greatest common divisor of N and 
yr/2 ±ij± quantum computer can find r efficiently by a se- 
ries of quantum operations on two quantum registers W and 
A. One is the work register W with L qubits, in which the job 
of finding the order is done; while the values of the function 
fix) are stored in the auxiliary register A with L' qubits. The 
sizes of the work and auxiliary registers are chosen as the in- 
tegers satisfying the inequalities TV 2 < q — 2 L < 2N 2 and 
2 L _1 < TV < 2 L . Here q is the Hilbert space dimension of 
the work register. 

As shown in figure 1, a realistic implementation of Shor's 
algorithm can be decomposed into the following unitary trans- 
formations: 

1) Initialize the work register in an equal-weight superpo- 
sition of all the logical states, and the auxiliary register in its 
logical ground state \0)a- Initially, each work qubit is in its 
logical ground state |0). Assuming that a Hadamard gate H 
is applied to each qubit in the work register at one time, the 
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|*(n)) = -= E e ~ lE]T1 b> ® e- lE ^ \0) A , (1) 

V 9 3=0 

with Ej being the energy of state \j) and h = 1. Here, 
T m (m — 1,2,3, ...) denotes the time interval between the 
(to — l)th and mth unitary operations. r+ = r m + e with 
e *C T m being the operational time of the mth unitary trans- 
formation, here assumed to be extremely small compared to 
other time scales. In other words, r+ refers to the time inter- 
val between the end of (m — l)th operation and the end of the 
mth operation. In what follows, the global dynamical phase 
exp(— iEqTi) will be omitted as it does not have any physical 
meaning. 

2) Calculate the function fN,a{j) = <F> mod N and then en- 
tangle the work and auxiliary registers \ f a ,N(s))A by 
applying a joint operation V. After another finite-time delay 
T2 before the next step (i.e., the third unitary transformation), 
the entangled state of the whole system becomes 

1 r_1 

|*(r 1 + +r 2 )) = — 53|V>w®|0>A, (2) 
V 1 s=0 

where 

\4>) A = CXp[-iE faJv(s) T 2 } \fa.N(s))A, 
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and 



3) Measure the auxiliary register \4>)a m its computational 
basis After this operation, the state of the whole sys- 
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V ex-p[-iE (lr+s) (T + + r 2 )l \lr + s) w , tem becomes l*( r i + T 2 )> = l^( T i + T 2 ))w ® l^i 
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with u> = [(5 — s — l)/r] being the largest integer less than 
(q — s — l)/r . The dynamical phases of the qubits in the work 
register, before and after the joint operation V, can be added 
directly, as this operator is diagonal in the logical basis. 
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T2 )}a- In other words, the work and auxiliary registers dis- 
entangle and the work register collapses to one of its periodic 



states \ip(i 
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For example, if the measurement on the auxiliary register 
1^)^ gives a value A s = a s mod N, then the work register 
immediately becomes 
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\^( T i + T i))w = — — = Vcxp[-i£: (; . r+s) (T 1 + + t+)] \lr + s) w - 
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After the third unitary transformation is applied, there is a third time delay 73. The state \tp(T^ + T^Yjw now evolves to 
I 



|V(T+ + T+ + T 3 )) W 
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exp [-~iE {lr+s) (T+ + t+ +T3)] \lr + s) w - 
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(3) 



Because of the collapse of the wavefunction ^(r-j*" + 72)) in 
Eq. (2), the dynamical phases accumulated by the wavefunc- 
tion \4>)a of the auxiliary register do not affect the algorithm 
anymore, as the relevant phase exp[— iEf a N r s \ t%] becomes 
a global phase. 

4) Perform the fourth unitary transformation: the quantum 
Fourier transform (^-Transformation) on the work register 
\ip)wt so that information regarding the order r of a mod N 
(i.e., the smallest integer r such that a r mod N = 1) can be 
more easily extracted. After the F-Transformation the state 
of the work register becomes 

\tP(t)) w =— Y,9(k)\k) w , 

V ? fe = 

with 

T = T++T++T+, 

being the time after applying the fourth unitary transforma- 
tion, and 
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After another delay time 74, i.e., right before applying the fifth 
unitary transformation, the work register evolves into 
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+ n)) w = 4= E 9(k) e- lEkTi \k) w . (4) 



5) Finally, we carry out a measurement on the work regis- 
ter in the computational basis {|j')vk} an d derive the desired 
order r satisfying the condition a r mod iV = 1. This mea- 
surement yields the state \k)w with probability 
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which is independent of the free evolution during the last de- 
lay 7-4. Notice that here P(k) only depends on the total effec- 
tive delay time r = t± + + 73 , but not directly on the 
individual time intervals r m , m = 1,2,3, 4. 

In this decomposition of Shor's algorithm we have included 
time delays only in between the various unitary operations, 
which were implemented by independently using various one- 
time evolutions ITMfl6ll . Note that only the delays from the 
initial Hadamard gates to the finishing Fourier transformation 
may result in physical effects. Fortunately, all the operators 
during these delays are either diagonal or at least not affecting 
the phase accumulation. Therefore, the phases in each qubit 
simply add up. 

If each unitary transformation is itself composed of several 
consecutive steps, with delays between these internal steps, 
we assume these delays to be negligible. This condition im- 
plies that the internal time delays occurring between steps 
within each unitary operation should be so short that their ac- 
cumulated phases are negligible. Such a condition is possibly 
difficult to satisfy experimentally. However, our results below 
show that even under such a restrictive condition the interfer- 
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ence effects due to dynamical phases between successive uni- 
tary transformation are already too significant to be ignored. 

For the ideal situation without any delay (r m = 0), the 
probability distribution P(k) in Eq. (5) reduces to that in the 
original Shor's algorithm |2]. However, Eq. (5) clearly shows 
that the expected probabilistic distribution may be strongly 
modified by the interferences due to the dynamical phases 
of the superposition wavefunction, which would consequently 
lead to a lower probability for obtaining the desired final out- 
put. 

III. EFFECTS OF DYNAMICAL PHASES 

In order to study the effects of dynamical phases, we need to 
compute P(k). The probability P(k) in (5) can be computed 
if the energies E^ r+s - ) for the various states \lr+s)w involved 
are known exactly. These will be computed below. 



A. Phase-matching condition for eliminating the coherent 
errors due to operational delays 



As a first approximation we assume that all qubits in a quan- 
tum computer system possess identical energy spectra. Such 
an approximation is valid for naturally identical systems like 
trapped ions. In this ideal case, when all the qubits have the 
same energy splitting between ground and excited states, dif- 
ferent quantum states with the same number of excited qubits 
will acquire the same dynamical phase. For example, the 
four-qubit states II3O2O1O0) and IO3O2O1I0) would acquire 
the same dynamical phase exp(— i3eot — i£\t) during a de- 
lay time t. Here eo and ei are the energies of a single qubit 
corresponding to the ground state |0) and the excited state |1), 
respectively. Under this approximation, equation (0 can be 
rewritten as 
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where A is the qubit energy splitting and Dp (L^) is the 
number of qubits in the logical state |1) (|0)) for the number 
state \lr + s)w- Obviously, when the total effective delay time 
t (r = + T2 + t^~) satisfies the phase-matching condition 



tA = {e 1 -e )T = 2mr, n = l,2,3,..., (7) 

the above probability distribution P(k) reduces to that of an 
ideal computation process with rA = 0. This implies that the 
interference due to the fast evolution of the dynamical phases 
can be suppressed periodically so that the correct results are 
obtained at the delay points indicated in |7}. 

Physically, this phase matching condition is related to the 
transformation of wavefunction from the interaction to the 
Schrodinger pictures. Theoretical derivations (see, e.g., 0) 
for realizing quantum computation are usually in the inter- 
action picture, in which the Hamiltonian for the qubit free- 
evolution does not appear, and the oscillation of the super- 
posed wavefunction does not exist. More specifically, if a 
system Hamiltonian H can be written as a sum of a free oscil- 
lator part and an interaction part H — Hq + V, so that the 
time-dependent Schrodinger equation can be written as (in 
the so-called Schrodinger picture where operators are time- 
independent while states evolve with time) 



one can introduce the interaction picture wavefunction 

|^s(i)) = exp(— iHot/H)] 1 ^/ (t)), which satisfies 

where Vi = exp(iHot/H)V exp(— iHot/H). Now that Hq has 
been eliminated from the Schrodinger equation, it seems that 
dynamical phases due to the qubit free evolution would have 
no effect. However, at the end of a calculation, physical mea- 
surements have to be performed to read out the computational 
results, and these measurements are generally performed in 
the lab frame (the Schrodinger picture), in which the dynami- 
cal phases reappear. More specifically, the measurement of an 
observable O can be expressed as (<i/)s(t)\0\ips(t)) — 
(xpi(t)\ ^(iHot/h)dexp(-iH t/h)\ipi(t)) 
(ipi (t)\0 1 (t)\ipi (t)) . In other words, if we prefer calculating 
the expectation value of a time-independent operator, it has to 
be done in the Schrodinger picture. If |V>/(V)) = J2j a j \j ) 
is the desired final state, the Schrodinger picture final state 
would take the form 

\Mr)) = E «J e ~^^> = E «J e - i(i ^° ))TA |i> • 
j 3 

(8) 

Therefore, the phase-matching condition @ would render the 
phases exp[— i(Lj — L^)tA] = 1, so that it enforces 
the equivalence of interaction picture and Schrodinger picture 
states, which ensures that the coherent error arising from the 
free evolution during the delay can be effectively eliminated. 
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In what follows, we illustrate our discussion with a few in- 
stances of Shor's algorithm. 

B. An analytical example for factoring a small composite 
number 

Let us first consider the factorization of the smallest com- 
posite number 4, which uses a two-qubit work register, a two- 



qubit auxiliary register and a — 3. After going through the 
four steps of Shor's algorithm as discussed above, the final 
work register state (Eq. (4)) is 



\lp(T + Ti)) W 
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(9) 



with C = 1 + e- lrA , and £ = 1 - e - lrA . Here, \a k ) w 
refers to the logical states (with a = 0, 1) of the fcth (with 
k = 0,1) qubit in the work register. In other hands, \0)w = 
|0i0 o )w, \l)w = |0il o ), \2)w = II1O0), and \3) w = 
|lilo)w- 

To derive Eq. (9), the measurement on the auxiliary reg- 
ister is the projection Pa = 1 1) a (1| a- Measuring the work 
register in the computational basis, the state (9) collapses 
to the expected one: either jQ)-^, or \2}w, with probability 
p e = |C| 2 — [1 + cos(rA)]/4. This implies that the desired 
results (|£| 2 = 1/2) are obtained, only if the phase-matching 
condition (7) is satisfied. Equation (9) also shows that the dy- 
namical phase acquired by each qubit after the Fourier trans- 
form does not result in any measurable physical effect. 



C. Numerical examples for factoring a few integers 

To quantitatively evaluate the effects of the dynamical 
phases when running Shor's algorithm, we introduce two 
delay-dependent functions: p e (k e ) is used to quantify the 
probability of obtaining the correct result k e , and 
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Pe = y^Pe(fce) 



(10) 



is the probability of computing all the correct outputs. P e = 1 
for an ideal computation process and for practical quantum 
computers at the phase-matching time intervals consistent 
with Eq. Q. For other delays not satisfying Eq. Q wrong 
results (k 7^ k e ) can be obtained so that P e < 1. 

We now run the algorithm to factorize N = 21 with a = 5 
using 9 work qubits. Fig. 2 shows the various outputs and the 
corresponding probabilities for different delay times r: tA = 
0, 0.47T, 7r, 1.67T, and 2ir. It is seen from Fig. 2 that, when 
the phase-matching condition Q is satisfied, the computed 
results are identical to that of an ideal computation process 
with rA = 0. Note in Fig. 2 that the maximum value of 
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FIG. 2: The probability P(k) (see Eq. (6)) of observing values of k 
for different values of rA = (ei — eo)r — 0, 0.4-7T, n, l.dnr, and 
2ir, given TV = 21, q — 512, a — 5, and the expected order r — 6. 
Here, r is the total effective delay time between unitary operations. 
The correct outputs are obtained when the phase-matching condition 
tA — 2n (or the ideal case r A = 0) is satisfied. The probabilities 
of obtaining the correct outputs far from the phase-matching condi- 
tions are very low (see the second, third, and fourth panels. Note 
the different scales for the vertical axes). Indeed, as shown in the 
bottom three panels, many incorrect results are produced when the 
phase matching condition given by Eq. Q is not enforced. 



P(k) w 0.2 at the matching condition and P(k) < 0.02 away 
from it. 

We plot the delay-dependent P e in Fig. 3 for several ex- 
amples: factorizing N = 15, 21, and 33, with a = 13, 5, 
and 5, and when using 4, 9, and 11 work qubits, respectively. 
As is shown in Fig. 3, the correct results are always obtained 
at the phase-matching time intervals given by Eq. (7). For 
other delay cases, especially near the delay points satisfying 
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FIG. 3: The probability P e of obtaining the correct results versus 
At = (ei — €q)t for running Shor's factoring algorithm in the pres- 
ence of delays. The lines for r = 4, 6, 10 correspond to the cases 
where 4, 9, 11 work qubits, given q — 16, 512, 2048, are used to 
factorize TV = 15, 21, 33 with a — 13, 5, 5, respectively. Note 
that the expected outputs can be obtained at phase-matching points: 
At = 2n, 4-7T. 



the condition: r A = (ej — e )r = (2n — l)ir, the correct 
results cannot be obtained (for the case where the expected or- 
der is a power of two; see, e.g., the continuous line for r = 4 
in Fig. 3) or may be obtained with very low probabilities P e 
(for the cases where the order r cannot divide the given q ex- 
actly; see, e.g., the lines for r — 6, 10 in Fig. 3). Of course, 
the dynamical oscillations can also be suppressed by trivially 
setting up individual delays r m as r A m = 2nir. The key ob- 
servation here is that only the total delay time, instead of the 
duration for every delay, needs to be set up accurately to avoid 
the coherent dynamical phase error. 

Classically, higher precision is usually obtained by using 
more computational bits. However, this is not necessarily the 
case in practical quantum computation. Indeed, in the cur- 
rent example of Shor's algorithm, after taking into consid- 
eration the influence of the time delays between consecutive 
computational operations, the more qubits are used, the lower 
the computational efficiency is. This relationship is clearly 
demonstrated in Fig. 4, which shows that the probability of 
obtaining any one of the correct results decreases exponen- 
tially when increasing the number of qubits of the work reg- 
ister. Such a scenario is to be expected, since the number of 
possible outputs in the final measurement increases exponen- 
tially with the number of the work qubits, which makes the 
constructive interference in Eq. (5) for the probability P{k) 
harder to achieve if rA deviates from the phase-matching con- 
dition (7). At the exact points when (ei — eo)r = 2nir, the 
constructive interference of the superposition wave functions 
ensures that the computational accuracy is independent of the 
number of involved qubits. 



D. Effect of energy splitting inhomogeneity 

In the calculations up to now, we have assumed that all 
qubits possess an identical energy splitting A = t\ — erj. In re- 
ality, especially for the solid state quantum systems such as the 
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FIG. 4: The probability p e (k e ) of obtaining one of the correct re- 
sults versus the number L of work qubits used to run the quan- 
tum algorithm factorizing TV = 15 in the presence of a delay 
At — (ei — eo)r = 57r/3. The straight line shows that this prob- 
ability p e (k e ) decreases exponentially with the number L of qubits 
used. The points on the line show the probability of obtaining one of 
the correct outputs k e = (0, 4, 8, 12) for 4-qubits, (0, 8, 16, 24) for 
5-qubits, (0, 16, 32, 48) for 6-qubits, (0 , 32, 64, 96) for 7-qubits, 
and (0, 64, 128, 192) for 8-qubits cases, respectively. 
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FIG. 5: The probabilities P e (for factorizing TV = 15 using 8 work 
qubits) of obtaining the correct results for different phase-matching 
cases: (A)t = 2-7T, 4-7T, 6-7T, 87t, with a common Gaussian energy 
splitting fluctuation with c/(A) = 0.5%. Note that this probability 
P e is higher at the phase matching points with shorter total delay 
time t. 



Josephson junction qubits and quantum dot trapped spins, dif- 
ferent qubits may have slightly different energy splittings due 
to system inhomogeneity. The logical states with the same en- 
ergy in the "identical qubit" assumption (e.g., II3O2O1O0) and 
IO3O2O1I0)) may now have slightly different energies. A crit- 
ical question then is how robust the phase matching condition 
Q is for a system of multiple qubits with fluctuations in the 
qubit energy splittings. Here we provide quantitative answers 
to this important question by numerically simulating Shor's 
algorithm assuming a Gaussian distribution for the qubit en- 
ergy splittings. In other words, the energy splitting Aj of 
the jth qubit is chosen randomly according to the distribution 
function 
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FIG. 6: The probabilities P e (for factorizing N = 15 using 8 work 
qubits) of obtaining the correct results for different fluctuations of 
energy splittings: a/(A) = 0.01%, 0.3%, 0.7%, 1.1%, with a com- 
mon phase-matching point: (A)r = 2n. Note that the probability 
at the pase-matching point is still sufficiently high, even if the en- 
ergy splittings of the qubits exist with certain fluctuations around the 
average value (A). 



around an average value (A) and width a. Figure 5 shows 
that the probability of obtaining correct answers decreases as 
the total time delay r increases. Also, Fig. 6 shows the depen- 
dence of P e on the width of the qubit energy splitting distri- 
bution (7, with the delay condition set at (A)r = 2tt. As ex- 
pected, a quantum computer runs with higher efficiencies for 
shorter time delays r and for narrower distributions -P(Aj) of 
energy splittings. In essence, here we study an effect similar 
to inhomogeneous broadening, which is not a true dephasing 
effect. This is consistent with our focus in this paper on the 
coherent errors instead of the incoherent ones. 



IV. CONCLUSIONS AND DISCUSSIONS 

When a real quantum computer performs a computational 
task, there must be unavoidable time intervals between con- 
secutive unitary operations. During these delays, the wave- 
function of a system with non-zero free Hamiltonian would 
acquire relative dynamical phases, if the two states for each 
qubit have different energies. These dynamical phases lead to 
fast oscillations in the total wavefunction, and modify the de- 
sired quantum interference required by quantum algorithms, 
which in turn reduce the probability of obtaining correct com- 
putational results. 

Here we have studied the effects of the dynamical phases in 
running a quantum algorithm (more specifically, Shor's fac- 
toring algorithm). We point out that a phase-matching con- 
dition can potentially help allieviate the interference prob- 
lems caused by the dynamical phases, and this condition is 
closely related to establishing the equivalence between quan- 
tum states in the Schrodinger picture and the interaction pic- 
ture through a quantum computation process. In the presence 
of coherent phase errors, we have demonstrated that, the prob- 
ability of obtaining the correct answer decreases exponentially 
with increasing number of qubits of the work register. In ad- 
dition, Shor's algorithm fails for the worst case scenario of 



tA = (2n — l)ir if the expected order r is a power of two. We 
have further shown that the phase-matching condition studied 
here is quite robust in the presence of small fluctuations in 
the qubit energy splittings. Unlike the refocusing technique 
in NMR experiments [3], which deals with unwanted evolu- 
tions due to uncontrolled qubit interaction, we have shown 
here that by properly setting the total effective delay, the un- 
wanted oscillations of the superposed wavefunctions due to 
the free Hamiltonians of the bare qubits can be effectively 
suppressed, thus the desired output can be obtained without 
additional operations. This implies that the quantum comput- 
ing may be performed in an effective interaction picture, in 
which coherent errors arising from the free evolution of the 
bare qubits during the operational delay can be automatically 
avoided. 

We emphasize that the present simplified approach only 
treats the delays between two sequential functional operations 
and neglects those inside these transforms. In fact, each func- 
tional transform, which is actually equivalent to a multi-qubit 
gate, can be, in principle, implemented exactly by using only 
one-time evolution fl5l Hal . This "coarse-grained" one-step 
implementation implies that the evolutions relating to the var- 
ious parts of the total Hamiltonian have been well controlled. 
Therefore, the operational delays, relating only to the free 
evolution ruled by the free Hamiltonian of the bare physi- 
cal qubits, within each one of these larger functional building 
blocks are assumed to be zero. Also, the dynamical phases ac- 
quired by the superposed wavefunctions can be added up for 
the operational delays before and after each functional trans- 
formation. Therefore, the phase-matching condition (7) exists 
for the total delay. 

The present calculation is done assuming that Shor's algo- 
rithm is accomplished in 5 lumped steps. A simple analysis 
can prove that, even if using an actual elementary gate array 
model, e.g., shown in figure 7 (for implementing the initializa- 
tions by using the Hadamard gates and the quantum Fourier 
transformation), the proposed phase-matching conditions (in 
terms of the total delay time instead of individual delay times 
of each operational delay) for avoiding the coherent phase er- 
rors are still valid. The key is that, only two elementary non- 
diagonal operations (i.e., R x in Hadamard gates) are applied 
to each qubit in the work register (see Fig. 7). The qubit is in a 
product state before the first non-diagonal R x gate, while the 
delays after the second non-diagonal R x in the corresponding 
Hadamard gate do not affect the results of projective measure- 
ment (see, e.g., Eq. (9)). Therefore, the dynamical phases 
acquired in different effective operational delays accumulate 
even when the operational delays inside the functional steps 
are considered. 

In the present approach, we have assumed that every qubit 
in the work register has the same waiting time Tj for each 
effective operational delay. In practice, this assumption is not 
necessary. Indeed, in the elementary gate array model, the 
waiting times for different qubits would have been different. 
However, the phase-matching condition Eq. (7) needs only a 
slight modification in this case, so that it becomes a condition 
for each qubit Q: r fc A fc = 2n k ir,k = 1,2, ...;n k = 1,2,.. 
for each qubit. Here, A k and r k are the energy splitting and 
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total controllable effective delay of the fcth qubit in the work 




register, respectively. 

Finally, we emphasize that the results presented in this pa- 
per, through obtained using a simple model for the delays, 
clearly demonstrate the necessity of taking into consideration 
the dynamical phases of the qubits in implementing quantum 
algorithms. 



(a) Hadamard gate (b) Quantum Fourier transformation 

FIG. 7: Quantum circuits formed by the elementary single- and 
two-qubit logic gates for performing (a) Hadamard gate for one 
qubit, and (b) quantum Fourier transformation for three-qubit. Here, 
Si (I — 1, 2, ...) and pki (k — 0, 1, 2, ...) refer to the operational de- 
lays inside them, respectively. In the logical basis, the single-qubit 
gate R z = exp(i7rcr z /4) and the two-qubit controlled-phase gate 
R k = |00)(00| + |01)(01| + |10)(10| +exp(2i7r/2 fc )|ll>{ll| are 
diagonal, while the single-qubit R x = exp(i7rcr :E /4) is not. 
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